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Project Goal 


The goal of this project is development of an optimization algorithm for use with 
nu m e ri cal silicon solar cell models. By coupling an optimization algorithm with a solar cell 
model it is possible to .imultaneously vary design variables such as impurity concentrations, 
front junction depth, back junction depth, and cell thickness to maximize the predicted cell 
efficiency. An optimization algorithm has been developed and interfaced with the Solar Cell 
Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to 
solve the differential equations which, along with several relations from the physics of sem- 
iconductors, describe m at h e m atically the operation of a solar cell A major obstacle is that the 
numerical methods used in SCAP1D require a significant amount of computer time, and dur- 
ing an optimization the model is called iteratively until the design variables converge to the 
values associated with the m axi mum efficiency. This problem has been alleviated by design- 
ing an optimization code specifically for use with numerically intensive simulations, to reduce 
the number of times the efficiency has to be calculated to achieve convergence to the optimal 
solution. Adapting SCAP1D so that it could be called iteratively by the optimization code 
provided another menu of reducing the cpu time required to complete an optimization. 
Instead of calculating the entire I-V curve, as is usually done In SCAP1D, only the efficiency 
is c alc u lated (maximum power voltage and current) and the solution from previous calcula- 
tions are used to initiate the next solution. Optimizations have been run for a variety of sub- 
strate qualities and levels of front and back surface passivation. This was done to determine 
bow these variables affect the optimized efficiency and the values of the optimized design 
variables. The sensitivity of efficiency to each of the design variables was investigated by 
changing one variable and reoptimizing the others. Work is progressing to include variables 
associated with the design of an anti-reflection coating in the optimization. 


Problem Statement 


How can all the inputs to a numerical model 
of a silicon solar cell be simultaneously 
varied to obtain the "optimal" design. 
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Outline 


I Define the Optimization Problem. 

a) Objective, variables, and constraints 

II Solving the Optimizatii i Problem 

a) Outline of optimization algorithm 

b) Calculating efficiency using SCAP1D 

c) Adapting SCAP1D for an iterative 
environment 

III Results of Optimization 

IV Future Work 


Doping Concentration Versus Position 
Use of Cerf Model for Doping 
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Optimization Model 


MAX Efficiency ( D 0 , Xp, D B , X B , D L , X L ) 

14 £ D 0 6 20.6 

14 <, D B < 20.6 

14 £ D L < 20.6 

0.1 Z X F <> 10.0 

0.2 £ X B S 50.0 

10.0 £ X L << 300.0 

0.0 ^ d b ~~ d b 

0.1 5 X L -Xp-X B 

Dq » log[Front surface doping concentration] P atoms/cm 3 
D b = log[Bulk doping concentration] B atoms/cm 3 

D l * logfBack surface doping concentration] B atoms/cm 3 
X F » Front junction depth urn 

X B = Back junction depth pm 

X B » Cel) thickness pm 
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SCAP1D (Solar Cell Analysis Program in 1 Dimension) 

Poisson's equation V J V « -2 ( n - p + N D - N A ) 

Electron Continuity Equation VJ n * q ( R - G ) 

Hole Continuity Equation ^J p * -q { R - G ) 

Hole snd Electron Carrier Transport Equations 

J„ ■ “9 m “ V( V + Y -^ ) + q D a Vn 

9 

* -9 Up P V( V - (1 -t) ) j - 9 D,, Vp 

where: 

effective gap reduction ■ A, * £ AEg + 0 O + 8 p j 

Ax + 

effective asymmetry factor * y 

Several assumptions are required for the validity of these equations. 

(1) Complete ionization of the dopants. 

(2) The heavily doped regions are quasi-neutral and in low iniection. 

(3) The parameters A, and y are functions of the earner concentrations. 

(4) Ag and y do not change when the device is not in equilibrium. 


I-V Curves 
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Power Versus Voltage 



Base Input Parameters 


Constant Input Paramours 


Q(!usuiiQoa 
Temperature 
Doping Profile 

Shadowing (including reflection) 
Auger Recombination 
Band Gap Narrowing 


lOOmW/cm 3 (AM 1.5) 

28 degrees C 
eric 

7% 

wmdA— H 

Slootboom DegraiT model 


inputs Vanad Parammcaily 
From surface recombination velocity (Sfl 
Back surface recombination velocity ($*) 

Minority carrier lifetime (used same formulas as in last laogress report) 
is electron minority earner lifetime. 

XpQ is hole minority earner lifetime (tlways taken as one half f 
For .2 ohm -cm substrate the input - 2 me, I ms, .4 mi gives bulk 
minority earner lifetimes of 54, 30, and 13 micro seconds respectively. 
#4 (back surfaoe reflection, 1.0 or 0.0) 


inputs Varied Paramotricaify or Opamistd 
From junction depth (Xf) 

Back junction depth (Xg) 

Cell thicknees (X[) 

From surface doping co nc en tr at ion (X>q) 
Bulk doping ooncemraaon (D B ) 

Back surface doping concentration {DQ 
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Optimal Solution for Case 2 


Parameters Held Constant During Optimization 


Front surface recombination velocity (Sf) 

1,000.0 

cm/s 

Back surface recombination velocity (S b ) 

1,000.0 

cm/s 

Electron minority carrier lifetime 1 (x^) 

1.00 

ms 

Hole minority carrier lifetime 1 (t^ 

0.50 

ms 

Optimal Values of Decision Variables 


Front junction depth (X F ) 

0.10 

| xm (lower bound) 

Back junction depth (X fl ) 

0.20 

|i m (lower bound) 

Cell thickness QC[) 

280.1 


Front surface doping concentration (D~) 

Zl53xi0 19 

P atoms/cm 3 

Bulk doping concemratioo (D B ) 

1.989x10“ 

B atoms/cm 3 

Back surface doping concentration (DJ 

3^14x10“ 

B stoms/cm 3 


Cell Performance Parameters 


Efficiency 

21.871 

% 

Open circuit voltage (V^) 

668.2 

mV 

Short circuit current dens (/*) 

39.063 

mA/cm 2 

Maximum power voltage (V^) 

384.86 

mV 

Fill factor 

0.8379 


Collection efficiency 

99.23 

% 

Bulk leaitivfy 

0.74 

ohm-cm 

Sheet rer'stance layer 1 

988.0 

ohm/O 

layer 2 

25 .37 

ohm/O 


1 Values in lightly doped silicon. 


Efficiency Versus Front Junction Depth 


Case 2 (t • l ms, sf » 1,000 cmJt, tb * 1,000 cm/s! 
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Optima! Front Surface Doping for Fixed Front Junction Depth 


Cav? 2 [T - 1 ms. « 1,000 cm/s, sb - 1,000 am/*] 



Efficiency Versus B&c.. Junction Depth 


Case 2 [t - 1 ms, if- 1.000 cm/s. sb - 1,000 an/s] 



Optima! Back Surface Doping for Fixed Back Junction Depth 

Cait 2 [t- 1 ms, if* 1 .000 cm/s, sb - 1,000 cm/*] 

a*-* — i 
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Efficiency Versus Back Junction Depth 


Case 3 [t - 1 ms, sf - 1,000 cm/s, sb - cm/s] 



Optimal Back Surface Doping for Fixed Back Junction Depth 


Coe 3 tt • 1 ms, sf • J ,000 cm/s, sb*~ cm/s] 
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Efficiency Versus r 

at Different Levels of Surface Passivation 


$ 



x 

1 « Intrinsic electron minority carrier lifetime 
sf [sb] - front [back] surface recombination velocity 
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Summary 


(1) Significant computational savings can be realized by adapt- 
ing a solar cell model for use with an optimization algo- 
rithm. 

(2) Comparisons should only be made between different values 
of a design variable after the other variables are optimized. 

(3) An optimization algorithm provides a systematic method 
for comparing different levels of technology and/or fabrica- 
tion processes. 

(4) A model coupled with an optimization algorithm provides a 
very powerful tool for analyzing the system modeled. 


Future Work 


(1) Analyze results of optimization runs. 

(2) Investigate other models for the doping concentrations. 

(3) Include a term for lateral resistance in the objective func- 
tion. 

(4) Include the design of the anti-reflection coating in SCAP1D 
and possibly in the optimization. 
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